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Abstract 

We provide an exact analytic description of decelerating, stopping and re-accelerating optical 
solitons in atomic media. By virtue of this solution we describe in detail how spatially localized 
optical memory bits can be written down, read and moved along the atomic medium in a prescribed 
manner. Dynamical control over the solitons is realized via a background laser field whose intensity 
controls the velocity of the slow light in a similar way as in the linear theory of electromagnetically 
induced transparency (EIT). We solve the nonlinear model when the controlling field and the 
solitons interact in an inseparable nonlinear superposition process. This allows us to access results 
beyond the limits of the linear theory of EIT. 
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Today, significant theoretical and experimental efforts are made to better understand the 
physics of light propagation in atomic vapors and Bose-Einstein condensates whose inter- 
action with light is well described by the nonlinear A-model. This process is accompanied 
by striking effects such as electromagnetically induced transparency which permits the 
propagation of light through an otherwise opaque medium; slow-light phenomena when the 
pulse is slowed down to velocities of a few meters per second yf, and optical information 

n n 

storage when the pulse is completely stopped inside the dielectric |3J,|4J. Even though the 
linear approach to describing these effects is developed in detail, it has reached its validity 
limits, because modern experiments require more adequate nonlinear descriptions. However, 
the theory of the nonlinear regime is still very incomplete. 

We investigate the problem of control over optical memory bits created by the slow-light 
solitons. Our theoretical model corresponds to a gas of alkali atoms, whose interaction with 
light is approximated by the structure of levels of the A-type (see Ref. |l|). The two lower 
levels are the states |1) and |2), while |3) is the excited state. The medium is described 
by the 3x3 density matrix p in the interaction picture. In order to cancel the residual 
Doppler broadening, two optical beams are chosen to be co-propagating. The fields are 
described by the Rabi-frequencies fl a ,b- The field Q a corresponds to a~ polarization, while 
the second fib corresponds to a + polarization. Within the slowly varying amplitude and 
phase approximation (SVEPA), dynamics of the atom- field system is well described by the 
reduced Maxwell-Bloch equations 
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Here ( = [z — z )/c, r = t — (z — z )/c, A is the detuning of the resonance, and v is 
the coupling constant. The matrix Hj = — | (fi a |3)(l| + fi 6 |3)(2|) + h.c. represents the 
interaction Hamiltonian. 

The system of equations Eq. (0) is exactly solvable in the framework of the Inverse Scat- 
tering (IS) method ^,0,3,0). We define the initial and boundary conditions for the system 
of equations Eq.(^Q) as follows. We consider a semi-infinite ( > active medium with a 
pulse of light incident at the point ( = (initial condition). This means that the evolution 



is considered with respect to the space variable (, while the boundary conditions should 
be specified with respect to the variable r. The boundary conditions are defined by the 
asymptotic values of the density matrix p and the matrix Hj at r — > ±00. 

In our previous work jij] we proposed a method to realize quantum memory bits in 
BEC by slow-light optical solitons. There we reported a slow-light soliton solution on a 
constant background field Qq, which was created by the auxiliary laser on entrance into 
the medium (£ = 0). The group velocity of the slow-light soliton depends explicitly on the 
field Qq. For the time-dependent background field some approximate solutions based on the 
methods of scaled time are studied in In a simplifying approximation the 

velocity is v g ~ c^. This expression immediately suggests that the soliton can be stopped 
through switching off the controlling laser field p, LL2f • I n the present paper we provide a 
full analytical description of how one can perform controlled preparation, manipulation and 
readout of optical memory bits formed in atomic vapors and Bose-Einstein condensates by 
slow-light solitons. 

Specifically, we consider a gas of rubidium atoms. To make parameters dimensionless, 
we measure the time in units of optical pulse length t„ = lfis typical for the experiments 
on the slow-light phenomena in rubidium vapors 3j, [l3| . The retarded time r is measured 
in microseconds and the Rabi frequencies are normalized to MHz. The spatial coordinate 
is normalized to the spatial length of the pulse slowed down in the medium, i.e. l p = v g t p « 
c 2u^v ^ n ex P er i men ts with rubidium the controlling background field is of order of a few 
megahertz. We choose VLq = 3. This corresponds to the group velocity of several meters per 
second, depending on the density of the atoms. We take the group velocity to be 10~ 7 c, so 
the pulse spatial length is 30/xm, and ( is normalized to 10~ 13 s. Then, in the dimensionless 
units, the coupling constant u = = 4.5. 

We consider the following scenario for the dynamics (see FigU). To begin with, in a 
remote past we create in the medium a slow-light soliton, and assume it is propagating on 
the constant background Qq. In the next step, we slow down the soliton by switching off the 
laser source of the background field. We assume an exponential decay of the background 
field with a decay constant a, i.e. f2o e ~ ar - At a certain, say T\ = 4/a, the field becomes 
negligible. Therefore, we cut off the exponential tail and approximate it by zero. At this 
step the soliton is completely stopped. The position where the soliton stops, depends on 
the decay constant a and on the moment when we switch off the laser The process of 
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FIG. 1: Time-dependence of the controlling background field. We choose a = 4, and the delay 
interval is T — T\ = 3. The dimensionless units are defined in the text. 

stopping the soliton corresponds to writing a memory bit into the medium. The information 
borne by the soliton is stored in the form of spatially localized polarization. This formation 
can live a relatively long period of time in atomic vapours or BEC. At the moment of time 
T we restore the slow-light soliton by abruptly switching on the laser. The whole dynamics 
is divided into four time intervals {J^ =0 T>i = (— oo, 0]U[0, T 1 ]U(T 1 , T]U(T, oo). The time- 
dependence of the intensity of the background field at entrance into the medium is given in 
Figffl 

Before the soliton enters the medium, the atoms are assumed to be in the state 
Notice that this state is a dark-state for the controlling field, which means that the atoms 
do not interact with the field created by the auxiliary laser. We specify the initial state of 
the system as follows: 

n(°) = o, ^ 0) = cXr), hU = |i>- (2) 

This configuration corresponds to a typical experimental setup (see e.g. The 
function Q(t) is given in FigO The state Eq.(j2J satisfies the Maxwell- Bio ch equations 
Eqs.(|H). Using the methods of our previous works Ji| [l^, we construct the one-soliton 
solution corresponding to the background profile (see FigJTJ), viz. 
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with the atomic state p = {ipa^i^atl, where 
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ltU = 11> + 



|2> 



2|A- A\w(r,Xy ' 2|A-A 



^at|3)- 



Here, 



7 , ^Cx 

' + ^ Im A-A 



^— + Re(z(r, A)) + In y/l + \w(r, A)p 



(4) 



0. = eo-^Re^^ + Im^A)), 

where A is an arbitrary complex parameter. The functions w(t, A), z(t, A) are of piecewise 
form, specific to each time region T>j. 



TABLE I: Exact analytical solution 
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For clarity we organize elements of the solution corresponding to different time regions 
in Table |U We use an auxiliary function w = f2 /(A + VA 2 + ^o 2 ), the index of Bessel 
functions is defined as 7 = (a + i\)/(2a). The values Cj of the constant C for each time 
region T>i are specified in the rightmost column of the table, the moment of time T is chosen 
as in FigJTJ i.e. T = 4/a + 3. Notice that in the table u>2 = w\(oo, A) and Z2 = 24(00, A). 
Therefore the solution in the region T>2 is parameterized by the asymptotic values of the 
data for the region D\ corresponding to the absence of cut-off of the exponentially vanishing 
tail. The region T> 2 describes the phase when the slow-light soliton is stopped, the fields 
vanish, while the information borne by the soliton is stored in the medium in the form of 
spatially localized polarization. At the moment of time T the laser is instantly turned on 
again. The stored localized polarization then generates a moving slow-light soliton. This 
process is described by the solution in the region V 3 . Except for the point Ti, the functions 
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w, z are continuous in r. This ensures that the physical variables such as the wave- function 
and field amplitudes evolve continuously. 
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FIG. 2: Contour plot of the intensity of Q a . We choose A = — 4.1i and zero detuning, A = 0. The 
break-up area in between the two solitonic trails manifests the creation of a standing memory bit 
in the medium. 

We demonstrate typical dynamics of the intensity of the field in the channel a in Fig|21 
This contour plot shows that in the process of rapid deceleration the solitonic trail forms 
end sharply. It is interesting to note that the restored pulse, however, has the same char- 
acteristics, i.e. the width and group velocity, as the input signal existed in the medium 
before the stopping. The dynamics of the corresponding localized polarization is presented 
in Fig|3] Notice that in the presence of the soliton the population flip from level |1) to |2) 
in the center of the peak is almost complete. A small fraction of the total population is 
located in the upper level |3) and provides for some atom-field interaction. In |5[ we show 
that the population of level |3) is proportional to \ flo\ 2 . Therefore when the field vanishes 
the population of the second level reaches unity, while the destructive influence of relaxation 
becomes negligible. This behavior of the system points to a sound possibility for realizing 
optical memory. 

To conclude, with an exactly solvable example, we have demonstrated a possibility to 
manipulate memory bits by slow-light solitons. We have proposed a method to prepare, 
control and read out optical memory bits in atomic vapors and Bose-Einstein condensates 
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FIG. 3: Population of level |2). Here, A = -4.L and A = 0. The time interval 1 < t < 4 
corresponds to a standing localized polarization flip. 

in the regime of strong nonlinearity. 
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